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GOAL 

The  FORTRAN  program  PFW.FTN  (Partially-Filled  Waveguide) 
was  written  to  compute  the  complex  propagation  constant  of  the 
hybrid  modes  propagating  in  an  inhomogeneously-filled  waveguide 
with  any  number  of  lossy  substrates  and  non-perfectly  conducting 
walls. 


FORMULATION 


The  general  version  of  the  program  solves  for  the  hybrid  TE  and 
TM  to  x  modes  that  propagate  in  structures  as  shown  in  Figure  1. 
To  determine  the  possible  propagation  constants  kz,  the  following 
transcendental  equations  corresponding  to  these  hybrid  modes  are 
solved  for  LSM  ( TMX )  modes  as 


t@Tl  kgi  — 


and  for  LSE  ( TEX )  modes 


— cotkx\h\  ~ 
Pi 


-  —  tankx2h2 
e2 


- COtkx  2^2 

/*2 


(i) 


(2) 


1 


with 

kl  +  i^  +  kl 

The  hybrid  modes  are  denoted  LSEmn  and  LSMmn.  The  tran¬ 
scendental  equations  have  an  infinite  number  of  solutions  for  a  given 
mode  number  n.  The  index  m  denotes  the  order  of  these  solutions. 
For  b  >  a,  the  dominant  mode  is  the  LSMoi  mode  (m  =  0,  n  =  1) 
followed  by  m  =  1,  and  so  on.  For  b  <  a,  the  first  mode  is  an  LSE 
mode  with  n  =  0.  In  this  case  the  numbering  on  the  index  m  starts 
at  m  =  1,  then  m  =  2,  and  so  on.  The  numbering  is  chosen  to  be 
consistent  with  that  of  the  empty  waveguide  [1]. 

The  equations  (1),  (2)  are  derived  using  the  transverse  resonance 
method,  since  the  structure  at  resonance  reduces  to  a  transmission¬ 
line  structure  with  propagation  in  the  x-direction.  In  the  case  of 
perfectly  conducting  walls,  the  end  transmission-lines  are  shortcir- 
cuited  as  shown  in  Figure  2.  However  finite  conductivity  of  the 
metal  housing  has  been  implemented  in  the  PFW.FTN  where  the 
top  and  bottom  walls  can  have  a  conductivity  a,  which  in  turns  can 
be  described  by  a  load  at  the  end  of  the  transmission  line.  The  pro¬ 
gram  also  accounts  for  substrate  losses  with  the  option  of  inputting 
the  dielectric  and  magnetic  loss  tangents. 

e,  =  (1  -  j  tan  <!>,)  (5) 

Hi  =  Hri  (1  j  tan  qfi).  (6) 

This  program  allows  to  compute  the  complex  propagation  con¬ 
stant  of  any  propagating  mode  in  a  rectangular  waveguide  struc¬ 
ture,  as  a  function  of  frequency.  It  can  also  calculate  the  cut-off 
frequency  of  the  modes  by  choosing  a  frequency  increment  small 
enough  around  kz  =  0. 

DESCRIPTION 

The  geometry  of  the  structure  is  inputted  through  a  user-friendly 
menu  program.  The  output  is  written  in  a  file  and  on  the  screen, 


(3) 

U)2C2l-l2- 

(4) 

O 


Figure  2:  Transmission-line  Analogy 


and  gives  the  propagation  constant  of  the  modes  propagating  within 
the  frequency  range  of  interest. 

.1  menu 

The  menu  is  an  interactive  program  that  allows  for  on-screen  cre¬ 
ation  of  the  input  data  file.  Also  changes  can  be  made  on  the  screen 
by  entering  the  variable  name  to  be  changed  (e.g.  FSTR  to  change 
the  starting  frequency).  The  program  will  then  prompt  for  a  new 
value  of  the  variable.  Dimensions  can  be  entered  either  in  centime¬ 
ters  or  inches.  Conductor  losses  in  the  top  and  bottom  walls  may 
be  included,  in  which  case  the  conductivity  of  the  walls  should  be 
given.  The  maximum  mode  number  n  refers  to  the  discussion  in  the 
previous  section.  The  output  is  written  in  a  file  called  ’RESULT’. 

.2  program 

The  PFW.FTN  program  calculates  the  roots  of  (1)  and  (2)  for  the 
lossless  case.  When  losses  are  considered,  the  roots  of  the  tran¬ 
scendental  equations  are  found  using  Muller’s  method  with  defla¬ 
tion.  The  present  program  calls  for  the  following  subroutines  of  the 
IMSL  library  (November  1,  1984  release)  [2] 

•  ZANLYT  (  find  the  zeros  of  a  univariate  complex  function  using 
Muller’s  method) 

•  UGETIO 

•  USPKD 

•  UERTST. 

These  subroutines  have  to  be  bound  to  the  main  program  PFW.FTN 
for  the  program  to  run.  If  this  library  is  not  available,  a  lossless  ver¬ 
sion  of  the  program  may  be  requested. 

.3  limitations 

The  program  PFW.FTN  is  a  generalized  code  that  can  handle  any 
number  of  dielectric  layers.  However,  for  programming  purposes, 
the  arrays  have  been  dimensionned  to  a  maximum  of  20  layers. 


4 


VERIFICATION 


Using  the  approach  described  in  the  previous  section  a  computer 
program  was  developed  to  calculate  the  complex  propagation  con¬ 
stant.  The  validity  of  this  program  has  been  verified  in  the  case  of 
lossless  substrates.  Good  agreement  with  Yamashita’s  paper  [3]  is 
shown,  and  the  results  of  the  program  PFW.FTN  given  in  the  next 
section  may  be  compared  to  Figure  3. 
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Freoucncy  (GHz) 


Figure  3:  Propagation  Constant  as  a  function  of  frequency  (comparison  with 
Yamashita) 


6 


I/O  FILES  AND  CODE 


.1  screen  sample 
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THIS  PROGRAM  CALCULATES  THE  PROPAGATION  CONSTANT  OF  THE  HYBRID  MODES 
EXCITED  IN  PARTIALLY-FILLED  WAVEGUIDES 

T.  EMI  LIE  VAN  DEVENTER  —  AND  —  LINDA  P.  B.  KATEHI 

MARCH  11,  1990  VERSION 


Fortran  PAUSE 
Type  return  to  continue 
Enter  name  of  configuration  data  file; 
yamashita 

ENTER  START  FREQUENCY  (GHz )  : 

10  • 

ENTER  STOP  FREQUENCY  (GHz)  : 

25 

ENTER  INCREMENT  FREQUENCY  (GHz) : 

1 

Select  units: 

(1)  inches 

[2]  centimeters 


Enter  #  OF  LAYERS  IN  WAVEGUIDE: 

2 

Enter  permittivity  for  layer  1 

1 

Enter  permittivity  loss  tangent  for  layer  1 

0 

Enter  permeability  for  layer  1 

1 

Enter  permeability  loss  tangent  for  layer  1 

0 

Enter  height  of  layer  1 
1.143 

Enter  permittivity  for  layer  2 
8.875 

Enter  permittivity  loss  tangent  for  layer  2 

0 

Enter  permeability  for  layer  2 

1 

Enter  permeability  loss  tangent  for  layer  2 

0 

Enter  height  of  layer  2 
.127 

Enter  waveguide  dimension  along  x: 

1.27 

Enter  waveguide  dimension  along  y: 

1.27 

Enter  maximum  #  of  modes  to  find 
3 

DIMENSIONS  ARE  TO  BE  ENTERED  IN  CENTIMETERS 


Name  of  geometry  configuration  file  (NAME]  yamashita 


the  units  are  CM 

cavity  dimension  along  x  (A) 

cavity  dimension  along  y  (B| 


1.2700000 

1.2700000 


START  FREQUENCY  TO  RUN  (GHz)  ( FSTR] : 
STOP  FREQUENCY  TO  RUN  (GHz)  (FSTP): 
INCREMENT  FREQUENCY  (GHz)  ( FNCR] : 


10.00000 

25.00000 

1.000000 


NUMBER  OF  DIELECTRIC  LAYERS  ( NLY ]  ;  2 


RELATIVE  ELECTRIC  PARAMETERS  OF  LAYERS: 

dielectric  loss 


Layer  1 
Layer  2 


constant 
[EL1R]  : 
[EL2R|  : 


1.000000 

8.875000 


tangent 
( ELI  1 1  : 
(EL2I)  : 


permeabl 1  ity 

real  imaginary 
0.000000  [ EL1M]  :  1.000000  (M1M)  : 
0.000000  ( EL2M)  :  1.000000  ( M2M)  ; 


0.000000 

0.000000 


GEOMETRY  OF  DIELECTRIC  LAYERS: 


CM 


Layer  1  thickness;  [LIB] 

Layer  2  thickness;  (L2B] 


1.143000 

0.127000 


(M| 


PROGRAM  PARAMETERS 
max.  n  mode  number 
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PACE  2 


03/12/90  1:33  PM 


//ter*/ uaera /▼•nd# renter /repo rta_dir/*cx*en . out 

CONDUCTOR  LOSSES 

Include  lower  ground  losses  :  [GLJ  F 

Include  upper  ground  losses  :  [Ul|  F 

Conductivity  of  the  walls  :  [SIGJ  0 . 00000000E*-C0 

Enter  Variable  Name  (**]  or  <re.  *jrn> 


operating 

frequency 

mode 

excited 

kz 

propagation 

kz/kO 

constant 
alpha  (Np/m) 

11.00000 

LSM_ 

0 

1 

31.31126 

0.13591 

O.OOOCO 

12.00000 

LSE 

1 

0 

79.74963 

0.31731 

C  .00000 

12.00000 

LSM_ 

0 

1 

120.35992 

0 . 4  7890 

0.00000 

13.00000 

LSE 

1 

0 

135.77371 

0.49867 

o  .ocooo 

13.00000 

LSM_ 

0 

1 

174 .46198 

0.64077 

0.00000 

14.00000 

LSE 

1 

0 

177.91298 

0.60677 

0,00000 

14.00000 

LSM_ 

0 

1 

223.21684 

0.76127 

0  .occoo 

15. OCOOO 

LSE 

1 

0 

214 . 72264 

0.68348 

0. ocooo 

15.00000 

LSM_ 

0 

1 

272 .61243 

0.86775 

o.coooo 

16.00000 

LSE 

1 

0 

248.85355 

0.74262 

0.00000 

16.00000 

LSE 

1 

1 

27.13668 

0.08098 

0.00000 

16.00000 

LSM 

0 

1 

326.30457 

0 . 973" 4 

o  .occoo 

16.00000 

LSM_ 

1 

1 

126.81892 

0.37845 

0.00000 

16.99999 

LSE 

1 

0 

281 . 65442 

0.79106 

0.  ocooo 

16.99999 

LSE 

1 

1 

1 34 . 67567 

0.37025 

o.coooo 

16.99999 

LSM 

0 

1 

386 . 91238 

1 .08669 

0.00000 

16.99999 

LSM_ 

1 

1 

188.03099 

0.52811 

0.00000 

17.99999 

LSE 

1 

0 

314.07828 

r .83312 

0  .ococo 

17.99999 

LSE 

1 

1 

193.52902 

0.51335 

0.00000 

17.99999 

LSM 

0 

1 

455.72186 

1.20884 

0.00000 

17.99999 

LSM 

1 

1 

233.84915 

0.62030 

0.00000 

17.99999 

LSM_ 

0 

2 

155.26567 

0.41186 

0.00000 

18 . 99999 

LSE 

1 

0 

347.03436 

0.87209 

0.00000 

18.99999 

LSE 

1 

1 

243. 39511 

0.61165 

0.00000 

18 . 99999 

LSM 

0 

1 

532 . 32367 

1.33772 

0.00000 

18.99999 

LSM 

1 

1 

272. 13550 

0. 68387 

0.00000 

18 . 99999 

LSM_ 

0 

2 

315. 90094 

0.79385 

O.OUOOO 

19.99999 

LSE 

1 

0 

381 .62134 

0.91105 

O.COOOO 

19.99999 

LSE 

2 

0 

87.68021 

0.20932 

0.00000 

19.99999 

LSE 

1 

1 

290.59106 

0.69374 

0.20000 

19.99999 

LSM 

0 

1 

614.82782 

1.46779 

C. 00000 

19.99999 

LSM 

1 

1 

305. 95987 

0.73043 

0 . ^0000 

19.99999 

LSM_ 

0 

2 

440.95157 

1.05269 

0 .oouuu 

20 . 99999 

LSE 

1 

0 

419. 31534 

0.95337 

0.00000 

20 . 99999 

LSE 

2 

0 

200.72998 

0.45639 

0.00000 

20.99999 

LSE 

1 

1 

338.57596 

0.76980 

0.00000 

20 . 99993 

LSM 

0 

1 

700.75470 

1.59327 

0.00000 

20 . 99999 

LSM 

0 

2 

554.51080 

1.26076 

0.00000 

20.99999 

LSM^ 

0 

3 

39.09147 

0.08888 

0.00000 

21.99999 

LSE 

1 

0 

461.96701 

1.00260 

0.00000 

21 . 99999 

LSE 

2 

0 

272.47369 

0.59135 

0.00000 

21 . 99999 

LSE 

1 

1 

390.15616 

0.04675 

0.00000 

21 . 99999 

LSE 

2 

1 

114.23763 

0.24793 

0.00000 

21 . 99999 

LSM 

0 

1 

787.95502 

1.71010 

0.00000 

21 . 99999 

LSM 

1 

1 

365.93430 

0.79419 

o.ooooc 

21 . 99999 

LSM 

2 

1 

71.16162 

0.15444 

0.00000 

21 . 99999 

LSM- 

0 

2 

661.28510 

1.43518 

0.00000 

21. 99999 

LSM_ 

0 

3 

362.40820 

0.78653 

0.00000 

2 2.99999 

LSE 

1 

0 

511.25302 

1.06133 

0.00000 

22. 99999 

LSE 

2 

0 

328.50262 

0.68195 

0.00000 

22 . 99999 

LSE 

1 

1 

447.42371 

0.92882 

0.00000 

22.99999 

LSE 

2 

1 

216.15340 

0.44872 

0.00000 

22. 99999 

LSE 

1 

2 

128.89127 

0.26757 

0.00000 

22 . 99999 

LSM- 

0 

1 

874 . 98932 

1.81642 

0.00000 

22.99999 

LSM 

1 

1 

393.56696 

0.81702 

0.00000 

22.99999 

LSM 

2 

1 

180.87822 

0.37549 

0.00000 

22.99999 

LSM- 

0 

2 

762.90967 

1.58375 

0.00000 

22.99999 

LSM_ 

0 

3 

525.42639 

1.09075 

0.00000 

23.99999 

LSE 

1 

0 

567.67902 

1 .12936 

0.00000 

23.99999 

LSE 

2 

0 

374.12955 

0.74431 

0.00000 

23.99999 

LSE 

1 

1 

510.94788 

1.01650 

0.00000 

23.99999 

LSE 

2 

1 

280.68002 

0.55840 

0.00000 

23.99999 

LSE 

1 

2 

278 . 37512 

0.55381 

0.00000 

23.99999 

LSM- 

0 

1 

961.05621 

1.91196 

0.00000 

23.99999 

LSM 

1 

1 

420.13651 

0.83584 

1 .78549 

23.99999 

LSM- 

2 

1 

244.50510 

0.48643 

0.00000 

23.99999 

LSM- 

0 

2 

060.26398 

1.71144 

0.00000 

23.99999 

LSM_ 

0 

3 

658.85938 

1.31076 

0.00000 

24 . 99999 

LSE 

1 

0 

630.29608 

1.20370 

0.00000 

24 . 99999 

LSE 

2 

0 

412.43860 

0.78770 

0.00000 

24.99999 

LSE 

1 

1 

579.72522 

1.10719 

0.00000 

24.99999 

LSE 

2 

1 

330.02109 

0.63029 

0.00000 

24 . 99999 

LSE 

1 

2 

390.52069 

0.74584 

0.00000 

24.99999 

LSM 

0 

i 

1045.78235 

1.99730 

0.00000 

24.99999 

LSM 

1 

1 

446.21091 

0.85220 

1.10039 

24 . 99999 

LSM- 

2 

1 

294  .  31281 

0.56210 

0.00000 

24 .99999 

LSM 

0 

2 

953.98407 

1.82198 

0.00000 

24 , 99999 

LSM 

1 

2 

124.01501 

0.23685 

0.00000 

24 . 99999 

LSM- 

0 

3 

777.25629 

1.48445 

0.00000 

Fortran  STOP 


FORTRAN  code 
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C 

c  t.e.  van  deventer 

c  radiation  laboratory 

o  room  3121,  eecs  bldg 

c  (313)  769-2975 

c  department  of  electrical  engineering  and  computer  science 

c  the  university  of  michigan 

c  ann  arbor,  mi  48105 

£****★*********************★*********★★★**********★★*★********★*★*★***★* 
c  program  pfw.ftn 

c  this  program  is  divided  into  3  parts, 
c  (1)  main  program 
c  (2)  subroutines 
c  (3)  functions 

c  date  :  o  31  may  1989 

£ 

c  bind  pfw.bin  muller.bin  -b  p 
c 

c  topic  : 

c  computes  the  propagation  constant  of  the  lse  and  Ism  modes 

c  to  solve  for  the  transcendental  equations  of  a  general 

c  n-layers  inhomogeneously-f illed  waveguide 

C 

£*********************************************************************** 

c - 

c  define  parameters 

c - 

integer  q, infer(l) ,n, nguess, nsig, kn, ms, ligne, nl 
integer  itmax, ier, i, 1, if ile, lunit, if,mflag 
integer  groundup, groundlow, miter 
integer  le,m, rdim,mid, r,mdim, index, ml, indexl 

real  a, b, h (20 ) , sig, ermax, mrmax, unit 

real  muO, epsO, pi, omega, kzr,rs(0:20),rs0(0:20),is(0:20) 

real  f op, fopstop, inert r, fopn 

real  epsr (20) ,mur (20) 

real  loss_tan (20) ,mu_tan (20) 

real  old_losstan (20 ) , old_mutan (20) 

real  inerkz 

real  rwi,  iwi, old_incrkz, eps, rkz,  fop_l 
real  attdbm, old_kzstop, as, old_incrf r 
real  absdet, abadetl, amp, differ 

complex  sum, sum_l, f sum 
complex  kz, kz_l, j, kO, kznoloss 
complex  xx ( 1 ) , zs, kzstart, cld_sum, eqn_c 
complex  er_c (20) ,mur_c (20) 

character*50  config_file 
character *2  hmode 

LOGICAL  do_gcl - do_ucl 
c - 
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c  common  blocks 

c - 

common/ cons t ant /muO, epsO , pi,  omega,  j ,  kO 

common/mike/er_c,  mur_c 

common /geomet ry/ zs, groundup, groundlow 

common/trick/ f sum 

common /dimens ions /mid 

common/param/mdim, le 

common/geom/a,  b 

common/ f req/ fop, fopstop, incrfr 

common/ layers/ rdim, h, eps  , mur , loss_tan, mu_tan 

common/operate/do_gcl , do_ucl, sig 

common /mode /m 

common/units/lunit, unit 

external  eq:i_c 

c - 

c  read  in  input  file 

c 

c  a  is  the  x-direction  (or  vertical) 

c  b  is  the  y-direction  (or  horizontal) 

c  h  is  the  height  of  the  layers 
c  fop  is  the  frequency  of  interest 
c  f op3top  is  the  maximum  frequency  of  interest 

c  incrfr  is  the  increment  of  the  frequency 

c  groundup  is  a  flag  (  1  -  losses  in  upper  wall  considered,  0  =  no  conductor  losses  ) 

c  groundlow  is  a  flag  (  1  -  losses  in  lower  wall  considered,  0  -  no  conductor  losses  ) 

c  incrkz  is  the  increment  of  the  propagation  constant  kz 
c  kzstop  is  an  upper  limit  to  the  value  of  the  propagation  constant 
c  if  le  is  1,  calculate  the  lse  propagation  constant 

c  if  le  is  2,  calculate  the  Ism  propagation  constant 

c 

c - 

call  logo 
call  menuc 

c - 

c  open  database 

c - 

open (unit-12,  file-' result' , status-'  unknown' ) 
call  writeout(12) 


groundup  -  0 
groundlow  -  0 
if  (do_ucl)  then 
groundup  -  1 
endif 

if  (do_gcl)  then 
groundlow  -  1 
endif 

c - 

c  change  values  to  SI  units 

c - 

fop  -  fop  *le9 
fopstop  -  fopstop  *  le9 
incrfr  -  incrfr  *  le9 


if  (lunit  .eq.  1)  then 


//tera/u3ers/vandeventer/moment_dir/pfw.ftn  Page3  03/12/90  12:13 

unit  =  .0254 

else  if  (lunit  .eq.  2)  then 
unit  “  le-2 
endif 

a  =  a  *  unit 
b  =  b  *  unit 
do  103  r=l,rdim 

h(r)  =  h(r)  *  unit 
103  enddo 

c - 

c  compute  some  constants 

c - 

pi  =  4 .  *  atan ( 1 . ) 
muO  =  4.e-7  *  pi 
epsO  =  le-9  /  (36.*pi) 
j  =  cmplx (0 . , 1 . ) 
mid  =  nint(rdim  /  2.) 

mrmax  =  0.0 
ermax  -  0.0 

do  101  r*=l,rdim 

er_c(r)  =*  epsr  (r)  *cmplx  (1 . 00,  -  loss_tan(r)) 
mur_c(r)  -  mur ( r) *cmplx ( 1 . 00, -  mu_tan(r) ) 
mrmax  -  amaxl (mur ( r)  ,  mrmax) 
ermax  =  amaxl (epsr (r) , ermax) 

101  continue 

ligne  -  0 

fop  -  fop  -  incrfr 

write  (12,*)  'operating  mode  propag' , 

&  'ation  constant  ' 

write  (12,*)  'frequency  excited  kz  kz/k0' 

s  , '  alpha  (Np/m) ' 

write  (12,  *)  ' - ' 

5  ,  ' - ' 

write  (*,*)  'operating  mode  propag', 

6  'ation  constant  ' 

write  (*,*)  'frequency  excited  kz  kz/k0' 

S  , '  alpha  (Np/m) ' 

write  ( *,  *)  ' - ' 

4  ,' - ' 

c - 

c  frequency  loop 

c - 

10  do  5  if-1,1000 

fop  -  fop  +  incrfr 

if  (fop  .gt.  fopstop)  then 
goto  5000 
endif 

fopn  -  fop  *  le-9 


omega  -  2.  *  pi  *  fop 
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kO  =  cmplx (omega  *  sqrt (muO  *  eps0),0.) 
kz  =  cmplx (. 01, 0 . )  *  kO 
incrkz  =  .001  *  kO 

miter  =  10000  *  (sqrt (ermax*mrmax) + . 2) 

if  (do_ucl  .or.  do_gcl)  then 
amp  =  sqrt (pi*fop*mu0/sig) 
else 

amp  =  0.0 
endif 

zs  =  cmplx (amp, amp) 

kzstart  =  kz 
old_incrkz  =  incrkz 

c - 

c  LSE  /  LSM  loop 

c - 

do  6  le-1,2 

if  (le  .eq.  1)  then 
hmode  «  ' E ' 

else  if  (le  .eq.  2)  then 
hmode  =  '  M' 
endif 

if  (le  .eq.  1)  then 
ml  -  0 
nl  -  1 

else 

ml  -  1 
nl  -  0 
endif 

c - 

c  m  -  mode  loop 

c - 

do  7  m-ml , mdim 

indexl  -  nl  -  1 
kz  -  kzstart 

17  continue 

incrkz  -  old_incrkz 
kz_l  -  cmplx (0 . 0, 0 . 0) 
sum_l  -  cmplx (0 . 0,  0 . 0) 
mflag  -  0 

c - 

c  loop  to  compute  the  kz  root  for  lossless  case 

c - 

12  do  20  ms-l, miter 

kz  -  kz  +  incrkz 

rkz  -  real(kz)  /  real(kO) 

if  (rkz  .gt.  sqrt (ermax'mrmax) * . 2 )  then 
c  write  <*,*)  'beta  >  quasi-static  value' 
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goto  1000 
endif 

sum  =  cmplx(0.,0.) 
call  hybrid (kz, sum) 


c  write  (*,*)  '  fop,m,  lekz, rkz, sum' , fopn,m, le, kz, rkz, sum 

c - refinement  loop - 

c 

c  the  loop  will  be  used  if  at  least  one  of  the  real  or  imaginary  parts 

c  of  the  sum  changed  sign  since  the  previous  case 

c - 

if  ( (real (sum_l)  .It.  0.0  .and.  real(sum)  .ge.  0.0) 

&  .or.  (real(sum_l)  .gt.  0.0  .and.  real(sum)  . le .  0.0) 

&  .or.  (aimag (sum)  .gt . 0 . 0  .and.  aimag (sum_l)  . It .  0 . 0) 

&  .or.  (aimag (sum) . It . 0 . 0  .and.  aimag (sum_l) .gt . 0 . 0) ) 

&  then 


mflag  »  mflag  +  1 
incrkz  *  incrkz  *  0.1 
kz  =  kz_l 

c - 

c  check  if  the  change  of  sign  is  due  to  a  singularity  in 

c  the  function  (mflag  test)  or  a  root  (differ) 

c - 

absdet  »  sqrt (real (sum) **2  +  aimag (sum) * *2) 
absdetl  -  sqrt (real (sum_l) **2  +  aimag (sum_l) **2) 
differ  -  abs (absdet  -  absdetl) 

if  (differ  .It.  10.)  then 
goto  13 
endif 

if  (mflag  .gt.  4)  then 
goto  17 
endif 

goto  12 
endif 

27  sum_l  -  sum 

kz_l  -  kz 

continue 
continue 
kznoloss  -  kz 
indexl  -  indexl  +  1 
rs(indexl)  ■  real(kz) 
rsO (indexl)  -  rkz 
is (indexl)  -  0.0 
c - 

c  loop  to  compute  the  kz  root  for  lossy  case 


c - initial  guess- 

40  rwi  -  real(kz) 

iwi  -  0.0 


xx(l)  •  cmplx (rwi, iwi ) 


20 

13 
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c - muller '  3  method  (determinant) 

50  continue 


eps  =  le-3 
nsig  =  3 
kn  =  0 
nguess  =  1 
n  =  1 

itmax  =  100 

call  zanlyt (eqn_c, eps, nsig, kn,  nguess, n,  xx,  itmax,  infer 
&  , ier) 

attdbm  =  -  aimag(xx (1) )  *  8.686 

rs(indexl)  -  real(xx(l)) 
rsO(indexl)  =  real(xx(l))  /  real(kO) 
is(indexl)  -  attdbm 

kz  =  kznoloss 
goto  17 

1000  continue 

do  153  ii=nl,indexl 
i  *  indexl+nl-ii 
ligne  =»  1 

write  (*,655)  fopn,hmode, ii,m, rs (i) , rsO (i) , is (i) 
write  (12,  655)  fopn,  hmode,  ii,m, rs (i) , rsO (i) , is (i) 

655  format (flO .5, 7x, ' LS' , Al, , i2, i2, 6x, fl2 .5, fl2 .5, fl2 .5) 

rs (i)  -  0.0 
rsO (i)  -  0.0 

i3  (i)  »  0 . 0 

153  enddo 

7  continue 
6  continue 

if  (ligne  .eq.  1)  then 

write  (12,  *)  ' - ' 

write  (*,*)  ' - ' 

endif 
ligne  •  0 
5  continue 

5000  continue 
stop 
end 

CS33333333333333333333333S3333333333S3333333S3333S3SSSSS3S3SSS333S3SSSS 
c  subroutines 

C3SS3333333333333333333333333333333333333333333333S33333S33333S33333333 
subroutine  hybrid (kz, sum) 


c  computes  the  arrays  corresponding  to  lse  and  Ism  common  terms 
c  le  -  (1)  corresponds  to  lse  modes  (electric) 


//tera/usars/vandaventar/momant_dir/pfw. ftn  Page  7 


03/12/90  12:13 


c  le  -  (2)  corresponds  to  Ism  modes  (electric) 

c 

c  - cannot  have  eps  or  mu  equal  to  zero  ! ! ! !- 

c 

c - 


integer  m,  le,  r,  rdim,mid,mdim 
integer  groundup, groundlow 

real  a,b 
real  ky 

real  h (20)  ,  epsr (20) ,mur(20) , loss_tan (20) ,mu_tan (2  0) 
real  muO , epsO ,  pi,  omega 

complex  k (20) , kx(20) , kz 
complex  j,k0,sum 
complex  er_c (20) ,mur_c (20) 
complex  zc (2, 20) , zO (2, 0 :2l) 
complex  twg, zs 

external  twg 

c - 

common/constant/muO, eps 0, pi, omega,  j ,  kO 
common/param/mdim, le 
common/mode/m 
c  ommo  n/mike/e  r_c , mu  r_c 
common/geometry/ zs, groundup, groundlow 
common/geom/a, b 

common/ layers /rdim, h, epsr,mur, loss_tan,  mu_tan 
c ommon / dimen s i on s /mi d 


ky  -  m  *  pi  /  b 
do  410  r*l,rdim 

k(r)  =  kO  *  csqrt (  er_c(r)*mur_c(r)  ) 
kx ( r)  »  csqrt (  k(r)**2  -  ky**2  -  kz**2) 
zc(l,r)  -  cmplx(omega  *  muO  *  mur_c(r))  /  kx(r) 
zc(2,r)  -  kx(r)  /  cmplx (omega  *  epsO  *  er_c(r)) 

410  continue 


i  •  le 

z0(i,0)  **  zs  *  groundup 

zO (i, rdim+1)  -  -zs  *  groundlow 

do  420  r*l,mid 

z0(i,r)  -  zc  (i,  r)  * 

&  (zO  (i,  r-1)  +  j*zc  (i,  r) *twg(kx (r) *h (r) ) ) / 

&  (zc(i,r)+j*z0(i,r-l) *twg (kx ( r ) *h (r) ) ) 

420  continue 


do  430  r«rdim, mid+1 , -1 
zO  ( i,  r)  -  zc  (i,  r)  * 

&  (zO(i,r+l)+j*zc(i,r)*twg(-kx(r)*h(r)))/ 
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&  (zc(i,r)+j*zO(i,r+l)*twg(-kx(r)*h(r))) 

430  continue 


sum  =  z0(i,mid)  -  z0(i,mid+l) 
c  write  (*,*)  '  kz,  i,mid,  zO (i,mid) , zO (i,mid+l) , sum' 

c  write  (*,*)  kz, i,mid, zO (i,mid)  ,  zO (i,mid+l)  ,  sum 

continue 

c - 

return 

end 

cssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss 
SUBROUTINE  MENUC 

INTEGER  i,j,il,i2,n 
INTEGER  menu 

INTEGER  rdim, mdim, le, lunit 

REAL  fop, fopstop, incrfr, sig, unit 

REAL  epsr (20) ,mur (20) , loss_tan (20) ,mu_tan (20) 

REAL  a, b 
REAL  x 
REAL  h (20) 

COMPLEX  eps (20) 

LOGICAL  do_gcl, do_ucl 
LOGICAL  exists 

CHARACTER* 80  warning 
CHARACTER* 10  check 
CHARACTER*8  ans,un 
CHARACTER* 50  config_file 
CHARACTER* 5 0  cfg_f ile, name_f ile 

CHARACTER*?  menu_elr (19) ,menu_eli (19) ,menu_elm(19) 

CHARACTER*?  menu_layer ( 1 9 ) ,menu_mu(19) 

********************************************************************************** 
common/param/mdim, le 
common/geom/a,  b 

common/f req/fop, fopstop, incrfr 
common/ layers/rdim, h, epsr,mur, loss_tan,mu_tan 
common/operate/do_gcl,do_ucl,  sig 
common/f iles/conf ig_f ile 
common/units/lunit, unit 

parameter (PI-3. 14159265358979324) 

data  menu_elr/' [EL1R] ' , ' [EL2R] ' , ' [EL3R] ' , ' (EL4R) ' , 

&  '  [EL5R] ' , '  [EL6R] ' , '  [EL7R] ' , '  [EL8R] ' , '  [EL9R]  '  , 

S.  '  [EL10R]  '  ,  '  [EL11R]  '  ,  '  [ EL12R]  '  ,  '  [EL13R]  '  ,  '  [EL14R]  '  , 

5  ' [EL15R] ' , ' [EL16R] ' , ' [ EL17R] ' , ' [EL18RJ ' , ' [EL19R] ' / 
data  menu_eli/' (ELI I] ' , ' [EL2I] ' , ' [EL3I] ' , ' [EL4I] ' , 

6  ' (EL5I ] ' , ' (EL6I ] ' , ' [EL7I] ' , ' [EL8I] ' , ' [EL9I] ' , 

5  ' [ELIO I ] ' , ' (ELI 1 I ] ' , ' (EL12I] ' , ' [EL13I] ' , ' [EL14I] ' , 

6  '  [EL15I] ' , '  (EL16I] ' , '  [EL17I] ' , '  (EL18I) ' , '  [EL19I]  '  / 
data  menu_elm/' [EL1M] ' , ' (EL2M] ' , ' [EL3M] ' , ' [EL4M] ' , 

S  ' [EL5M] ' , ' [EL6M] ' , ' [EL7M] ' , ' [EL8M] ' , ' [EL9M] ' , 
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&  ' [EL10M] ' , ' [EL11M] ' , ' [EL12M] ' , ' [EL13M] ' , ' [EL14M] ' , 

&  ' [EL15M] ' , ' [EL16M] ' , ' [EL17M] ' , ' [EL18M] ' , ' [EL19M] ' / 

data  menu_mu/ ' [M1M] ' , ' [M2M] ' , ' [M3M] ' , ' [M4M] ' , 

&  ' [M5M] ' , ' [M6M] ' , ' [M7M] ' , ' [M8M] ' , ' [M9M] ' , 

&  '  [M10M]  '  ,  '  [MUM]  '  ,  '  [M12M]  '  ,  '  [M13M]  '  ,  '  [M14M]  '  , 

&  '  [M15M] ' , '  [Ml 6M] ' , '  [M17M] ' , '  [M18M] ' , '  [M19M] ' / 

data  menu_layer/' [LIB] ' , ' [L2B] ' , ' [L3B] ' , ' [L4B] ' , 

&  '  [ L5B] ' ,  '  [L6B] ' , '  [L7B] ' ,  '  [L8B] ' , '  [L9B] '  , 

&  '  [L10B] ' ,  '  [L11B] ' ,  '  [L12B] ' ,  '  [L13B] ' , '  [L14B] ' , 

&  '  [L15B] ' ,  '  [L16B] ' ,  '  [L17B] ' ,  '  [L18B] '  ,  '  [L19B] '  / 


DEFINE  WRITES  A  FILE  CONTAINING  ALL  OF  THE  PARAMETERS  NEEDED  FOR 
THE  ANALYSIS  OF  THE  INHOMOGENEOUSLY-FILLED  WAVEGUIDE  PROBLEM. 

THE  VARIABLES  AND  PARAMETERS  ARE  NAMED  AS  FOLLOWS: 


a 

b 

fop 

fopstop 
incrf r 
rdim 


dimension  of  cavity  along  x  axis 
dimension  of  cavity  along  y  axis 
frequency  of  analysis  in  GHz 
maximum  frequency 
frequency  increment 

number  of  layers  in  cavity  substrate 


C  THE  NEXT  SEGMENT  DETERMINES  WHETHER  AN  EXISTING  CONFIGURATION  FILE  IS 

C  TO  BE  USED 

A************************************************,,******************************** 

100  write (*,*)  'Enter  name  of  configuration  data  file;  ' 

read (*, ' (A50) ' )  config_file  [geometric  configuration  file  name 

cfg_file=conf ig_f ile 

inquire (f ile=cfg_f ile, exist=exists)  (check  if  format  file  already  exist 
call  trim (cfg_f ile, 50, il, i2) 
if  (exists  )  then 

warning  =  '  WARNING:  File  ' / /cfg_f ile ( il : i2 ) / / 

&  '  already  exists.  Use  existing  file?'// 

&  '  (  <return>  =  yes  ) ' 

write  (6,*)  warning 

read  (5, '  (al) ' )  ana 

if  ( (ans . ne . ' y' ) . and. (ans . ne . ' Y' ) . and . (ans . ne . ' ' ) )  goto  100 
goto  200 
endif 


********************************************************************************** 
C  NEXT  SECTION  OF  CODE  IS  ONLY  REACHED  IF  NO  CONFIGURATION  FILE  ALREADY 

C  EXISTS  AND  WILL  COLLECT  ALL  OF  THE  NECESSARY  INPUT  AND  GENERATE  A  NEW 

C  CONFIGURATION 

********************************************************************************** 

949  write <*,*)  'ENTER  START  FREQUENCY  (GHz)  :  ' 
read(*, *,err-949)  fop 

950  write (*,*)  'ENTER  STOP  FREQUENCY  (GHz)  :  ' 
read ( *, *, err-949)  fopstop 

951  write ( * , * )  'ENTER  INCREMENT  FREQUENCY  (GHz):  ' 
read ( *, *, err«949)  incrfr 

write(*,*)  'Select  units:  ' 
write<*,*)  '[1]  inches  ' 
write(*,*)  '[2]  centimeters  ' 
read ( *, *)  lunit 
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write ( * , * ) 

write (*,*)  'Enter  #  OF  LAYERS  IN  WAVEGUIDE:  ' 
read(*,*)  rdim 

x  =  0.0 

do  i  =  1,  rdim 

write (*, *) 
write (*, *) 

&  'Enter  permittivity  for  layer  ',i 

read(*, *)  epsr (i) 

if  (epsr(i)  .eq.  0.0)  epsr(i)  =  1.0 
write (*, *) 

&  'Enter  permittivity  loss  tangent  for  layer  ',i 

read(*,*)  loss_tan(i) 
write (*, *) 

&  'Enter  permeability  for  layer  ',i 

read(*,*)  mur(i) 

if  (mur(i)  .eq.  0.0)  mur(i)  =  1.0 
write  (*, *) 

&  'Enter  permeability  loss  tangent  for  layer  ',i 

read(*,*)  mu  tan(i) 


write (*,*)  'Enter  height  of  layer  ',i 
read  (.*,*)  h(i) 

x  -  x  +  h  (i) 
enddo 

********************************************************************************** 

*  WAVEGUIDE  GEOMETRY 

********************************************************************************** 
write (*,*)  'Enter  waveguide  dimension  along  x:  ' 
read(*, *)  a 

write (*,*)  'Enter  waveguide  dimension  along  y:  ' 
read(*,  *)  b 

if  (x  ,ne.  a)  then 
goto  499 
endif 

*****************************************  ***************************************** 

*  NUMERICAL  PARAMETERS 

********************************************************************************** 

write (*,*)  'Enter  maximum  #  of  modes  to  find' 
read ( *, *)  mdim 

dogcl  -  .false. 
do_ucl  -  .false, 
le  -  2 
sig  -0.0 


goto  400 
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£★*****★****★★★★★*★★★★★★★****★★***★★★★**★*★★★★★★★★★********★★**•*★***★********** 

*  THE  NEXT  SEGMENT  READS  DATA  FROM  EXISTING  FILE 

q****************************************************************************** 

200  open (10, f ile=cfg_f ile) 
read  (10, *) 

read  (10,650)  config_file 

read  (10, *) 

read  (10,*)  mdim 

read  (10,*) 

read  (10,*)  lunit 

read  (10, *) 

read  (10, *)  a,b 

read  (10, *) 

read  (10,*)  f op, f opstop, incrf r 

read  (10, *) 

read  (10,*)  rdim 

read  (10,*) 

do  i  =  1 , rdim 

read (10, *)  h (i) , epsr (i) , mur ( i) , loss_tan ( i) , mu_tan (i) 
enddo 

read  (10,*) 

read  (10,*)  do_gcl 

read  (10,*)  do_ucl 

read  (10,*) 

read  (10,*)  sig 

read  (10,*) 

read  (10,*) 

650  format (A50) 

clo3e (10) 

C  THE  NEXT  PART  IS  THE  MAIN  PART  OF  THE  MENU  SYSTEM 
************************************************************************ ********** 

•I******************************************************************************** 

*  WRITE  VARIABLE  LIST  ON  SCREEN 
********************************************************************************** 

400  continue 

if  (lunit. eq. 2)  then 

write  (*,*)  'DIMENSIONS  ARE  TO  BE  ENTERED  IN  CENTIMETERS  ' 
un  -  'CM' 

else  if  (lunit. eq.l)  then 

write (*,*)  'DIMENSIONS  ARE  TO  BE  ENTERED  IN  INCHES  ' 
un  -  ' INCH' 
endif 


write)*, *) 

write ( *, *)  '************ 
write ( *, *) 

write (*,705)  config_file 
write (*, *) 
write)*, 704)  un 
write)*, 710)  a 
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write (*,  715)  b 
write  ( * ,  *) 


704  format ('the  units  are  ',A7) 

705  format ('Name  of  geometry  configuration  file 


[NAME] 


A50 ) 


710 

format (' cavity  dimension 

along  x 

[A] 

.  r 

•  t 

& 

f 14 . 7) 

715 

format (' cavity  dimension 

along  y 

IB] 

,  r 
•  t 

& 

f  1 4 . 7 ) 

write)*,*)  'START  FREQUENCY  TO  RUN 

(GHz) 

[FSTR] :  ' 

,  fop 

write)*,*)  'STOP  FREQUENCY  TO  RUN 

(GHz) 

[ FSTP ]  :  ' 

, fopstop 

write)*,*)  'INCREMENT 
write (*,  *) 

FREQUENCY 

(GHz) 

[ FNCR ] :  ' 

, incrfr 

716 

format ( ' frequency'  ,  i3,  ' 

in  GHz' , 3x, A7 

, 2x, f 1 4 . 7 ) 

write (*, *) 

write)*,*)  'NUMBER  OF 

DIELECTRIC  LAYERS 

[NLY]  :  ', 

rdim 

if  (rdim.gt.l)  then 

write (*, *) 

write)*,*)  'RELATIVE  ELECTRIC  PARAMETERS  OF  LAYERS:' 

write ( *, * )  ' 

dielectric 

f 

t 

& 

'  loss 

permeability' 

write ( * , * )  ' 

constant 

t 

t 

& 

'  tangent 

real 

/ 

t 

& 

'  imaginary' 
do  i  «  1 , rdim 

write ( *, 720 )  i, menu_elr ( i) , epsr ( i) , 

menu_eli (i) 

r 

& 

loss_tan  < i) 

,menu_elm)i)  , 

mu  r ( i ) 

, menu_mu ( i ) 

r 

& 

mu  tan ( i) 

enddo 


write ( *, *) 
write(*,7l8)  un 
write (*,  *)  ' 


do  i  -  1 , rdim 

write(*,725)  i, menu_layer ( i) , h ( i) 
enddo 


write  (*,  *) 
endif 

718  format ('GEOMETRY  OF  DIELECTRIC  LAYERS:  ',2x,A7) 

72  0  format ('Layer' , i3, 5x,A7, '  : ' , f 12 . 6, 2x, A7,  ' :',fl2.6, 

S  2x, A7 , ' :',fl2.6,2x,A7,' : ' , f 12 . 6) 

725  format (' Layer' , i3,  '  thickness;  ' ,  A7, '  : ' ,  f 12 . 6) 


write (*, *) 

write)*,*)  'PROGRAM  PARAMETERS  ' 
write)*, 777)  mdim 
write ( *, *) 

write ( * ,  *)  'CONDUCTOR  LOSSES  ' 
write)*, 744)  do_gcl 
write(*,745)  do_ucl 
write**, 746)  sig 


777  format ('max.  n  mode  number 

744  format (' Include  lower  ground  losses 


(M)  ',  i 5 ) 

8x, ' [GL]  ' , 11, 5x> 
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745  format (' Include  upper  ground  losses  :',8x,'[UL] 

746  format (' Conduct ivity  of  the  walls  :',8x,'[3IG) 

c  goto  9999 

★★fr******************************************************* 

***************★★**★***■*■★**★★**★★*****★*■****★********★*★** 

A-********************************************************* 

C  THE  NEXT  PART  READS  THE  RESPONSES  TO  THE  MENU  PROMPT 

C  DETECTS  WHICH  VARIABLE  IS  TO  BE  CHANGED 

C  PROMPTS  FOR  THE  CHANGE  AND  MAKES  IT 

C  AND  UPDATES  THE  MENU 

**•*★**★★***★*★***★★*********★*********•★**★*★*★**★*★★★*★★★★ 


write ( * , * ) 

write'*,*)  'Enter  Variable  Name  [**]  or  <return>' 
read{*,' (A4) ' , err=499)  check 
********  ALTER  VARIABLES  AS  DESIRED 


if  (check . eg )  goto  600 


if 


& 

4 

4 


<  (check .eq. ' STOP' ) 
. or .  ( check . eq . ' QUIT’  ) 
.or.  (check .eq. ' EXIT'  ) 
.or. (check . eq. ' ABORT' 


.or.  (check . eq .' stop' ) 
.or.  (check . eq .' quit ' ) 
.or.  (check . eq .' exit ' ) 
.or.  (check . eq. ' abort '  )  ) 


if  ( (check . eq. ' name' ). or . (check .eq .' NAME' ) )  then 
write  (*,*)  'ENTER  NEW  CONFIGURATION  FILE  NAME  ’ 
read(*, '  <a50) ' ,err»499)  name_file 
call  t r im (name_f ile, 50 , il ,  i2 ) 
inquire ( f ile-name_f ile (il : i2) , exist=exists ) 
if  (exists)  then 

warning  =  '  WARNING:  File  ' //name_f ile (il : i2) / 
4  '  already  exists.  Input  existing  file? 

4  '  (  <return>  =  yes  ) ' 

write  (6,*)  warning 
read  (5, '  (al) ' )  ans 

if  ((ans.ne.'y')  .and.  (ans.ne.'Y')  .and.  (ans.ne." 
conf ig_f i le»name_f ile 
cfg_f ile»name_f ile 
goto  200 
el  se 

C3"f ig_f ile-name_f ile 
cfg_f ile«name_f ile 
goto  400 
endif 
endi  f 


if  ( ( check . eq. ' a ' ) . or . (check . eq . ' A' ) )  then 
write ( », *)  ' ENTER  A  ' 

read ( *, *, err-499)  a 
goto  400 
endif 

if  ( (check  eq .’ b' ). or . (check . eq .' B' ) )  then 
write  (  *, »)  'ENTER  B  ' 


' , 11, 5x) 

'  ,  el5 . 8) 

★  *****★★*★*★*★★★★★*•*•**** 
A*********************** 
************************ 


**************  **■*★★**•*** 


stop 


/ 

'  // 


) )  goto  501 
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read ( * , * , err=499)  b 
goto  40'' 
endif 

if  ( (check . eq f st r' ). or . (check . eq FSTR' ) )  then 
write (*,*>  'ENTER  START  FREQUENCY  (in  GHz)' 
read (*,*, err-499)  fop 
goto  400 
endi  f 

if  ( (check . eq .' fstp' ). or . (check . eq .' FSTP ') )  then 
write  (*,*)  'ENTER  STOP  FREQUENCY  (in  GHz)' 
read ( * , * , er r=49 9 )  topstop 
goto  400 
endif 

if  ( (check . eq .' fncr' ). or . (check . eq. ' FNCR' ) )  then 
write  (*,*)  'ENTER  FREQUENCY  INCREMENT  (in  GHz)' 
read (*,*, err-499)  incrfr 
goto  400 
endif 

if  ( (check . eq .' nly' ). or .  (check . eq .' NLY' ) )  then 
write  (*,*)  'ENTER  #  r>F  DIELECTRIC  LAYERS' 
read ( " , ", err=499)  rdim 
goto  400 
endif 

do  i  -  1,  rdim 

if  ( (check ( 1 : 1) .eq .' e' ). or . (check (1 : 1) . eq. ' E' ) )  then 
if  < (check  (2 : 2 )  . eq . ' 1' )  . or .  (check (2 : 2)  . eq. ' L' ) )  then 

if  < (check ( 4 : 4)  . eq. ' r' )  . or .  (check ( 4 : 4 )  .eq . ' R' ) )  then 
if  (check  ( 3  :  3)  ,eq.menu_elr ( i)  ( 4  :  4) )  then 
write (*,*)  'ENTER  ' , menu_elr ( i ) 
read ( *, *, err-499)  epsr(i) 
eps(i)-  epsr(i)  *  cmplx ( 1 ., loss_tan (i ) ) 
goto  400 
endif 

else  if  (  (check  (4:4)  . eq . ' i ' )  .or.  (check(4:4)  . eq . ' I' ) ) 
if  (check (3 : 3)  .eq.menu_eli (i)  (4 : 4) )  then 
write  (*, *)  'ENTER  ' , menu_eli ( i) 
read (*,*» err-499)  loss_tan(i) 
eps(i)  -  epsr(i)  *  cmplx ( 1 ., loss_tan ( i) ) 
goto  400 
endif 

else  if  ( (check ( 4 : 4) .eq. ' m' ) . or . (check ( 4 : 4 ) . eq . ' M' ) ) 
if  (check (3 : 3)  .eq .menu_elm ( i)  ( 4  :  4) )  then 
wr4te(*,*)  'ENTER  ' ,menu_elm(i) 
read (*,*, err-499)  mur(i) 
goto  400 
endif 

else  if  ( (check (5 : 5) .eq. ' r' ) .or . (check (5 : 5) .eq. ' R' ) ) 
if  (check  (4 : 4)  ,eq.menu_elr (i)  (5  :  5) )  then 
write (*#*)  'ENTER  ' , menu_e 1 r ( i ) 
read(*, *, err-499)  epsr(i) 
eps(i)-  epsr(i)  *  cmplx ( 1 ., loss_tan ( i) ) 
goto  400 
endif 


then 


then 


then 
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else  if  ( (check (5:5)  .eq. ' i' )  .or  .  (check (5: 5)  .eq. ' I' )  )  then 
if  (check  (4  :  4)  .eq.menu_eli  (i)  (5  :  5) )  then 
write!*, *)  'ENTER  '  ,  menu_e li  (  i ) 
read ( * , * , err=4 99)  loss_tan(i) 
eps(i)  =  epsr(i)  *  cmplx ( 1 . , ioss_tan ( i ) ) 
goto  400 
endif 

else  if  ( (check  ( 5  :  5)  . eq . ' m' )  . or .  (check  (  5  :  5)  . eq . ' M' ) )  then 
if  (check<4:4) .eq. menu_elm ( i ) ( S : 5) )  then 
write  (*,*)  'ENTER  ' ,menu_elm(i) 
read ( * , * , e r r=4 99 )  mur(i) 
goto  400 
endif 
endif 
endif 
endif 

if  (  (check ( 1 : 1 ). eq .' 1' ). or . (check ( 1 : 1 ). eq .' L' )  )  then 

if  (  (check ( 3 : 3)  . eq. ' b' )  . or .  (check  ( 3  :  3)  .eq . ' B' )  )  then 

if  (  check  (2 :2)  .eq.menu_layer (i)  (3  :  3)  )  then 
write(*,*)  'ENTER  ' , menu_layer (i) 
read  ( * ,  * ,  er  r**4  99)  h(i) 
goto  400 
endif 

else  if  ( (check  ( 4  :  4)  . eq. ' b' )  . or .  (check ( 4 : 4 )  . eq. ' B' ) )  then 
if  (  check ( 3 : 3)  . eq . menu_layer ( i)  ( 4  :  4 )  )  then 

write!*, *)  'ENTER  ' , menu_layer ( i ) 
read (*,*, er r-499 )  h(i) 
goto  400 
endif 
endif 
endif 

if  (  (check  ( 1  : 1 )  . eq . ' m' )  . or .  (check  (1:1)  . eq . ' M' )  )  then 
if  (  (check ( 3 : 3)  . eq . ' m' )  . or .  (check ( 3  :  3)  .eq. ' M' )  )  then 
if  (  check  (2  :  2)  ,eq.menu_mu (i)  (3 : 3)  )  then 
write  (*,*)  'ENTER  ' ,menu_mu(i) 
read ( * , * , err-499)  mu_tan(i) 
goto  400 
endi  f 

else  if  ( (check ( 4 : 4 )  .eq. ' m' )  . or . (check ( 4 : 4 )  . eq. ' M' ) )  then 
if  (  check  ( 3  :  3)  . eq ,menu_mu ( i)  (  4  :  4 )  )  then 

write!*,*)  'ENTER  ' ,menu_rau(i) 
read(*, *, err-499)  mu_tan(i) 
goto  400 
endif 
endif 
endif 
enddo 


if  ( (check . eq .' m' ). or . (check .eq .' M' ) )  then 
write ( *, •)  'ENTER  M' 
read (*,*, er r-499)  mdim 
goto  400 
endif 

if  ( (check . eq .' gl' ). or . (check .eq .' GL' ) )  then 
do_gci  -  .not.do_gcl 
goto  400 
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endif 

if  ( (check .eq. '  ul' ). or .  (check . eq UL' ) )  then 
do_ucl  *  .not.do_ucl 
goto  400 
endif 

if  ( (check .eq. ' sig' ) .or . (check .eq. ' SIG' ) )  then 

write (*,*)  'ENTER  CONDUCTIVITY  OF  THE  WALLS  (S/m)' 
read ( * ,  * , err=4 99 )  sig 
goto  400 
endif 

499  write(*,*) 
write (*, *) 
write  (*,  *) 
write ( * , * ) 
write  ( * ,  * ) 
goto  400 


*********************************************************************#**#****#ft*** 

C  THE  NEXT  PART  IS  THE  LAUNCHING  POINT  FOR  THE  REST  OF  THE  PROGRAM 

C  IT  IS  ONLY  REACHED  WHEN  FINISHED  WITH  THE  MENU  PART 

ft********************************************************************************* 

********************************************************************************** 

600  continue 

********************************************************************************** 
********  WRITE  VARIABLES  TO  NEW  FILE  OR  OLD  FILE 
500  call  trim (cfg_f ile, 50, il, i2) 

inquire ( f ile-cfg_f ile, exist -exists) 
if (exists )  then 

write(*,*)  'File  already  exists  - 

&  'overwrites  not  allowed.' 

write  (*,*)  'Rename  [R]  or  <return>  ( —  changes  not  saved)' 
read ( *,  ' (A8) ' )  ans 
if  (ans.ne."  )  then 

if  (ans . ne .' override' )  then 
goto  501 
else 

goto  599 
endif 
endif 
else 

599  open  ( 1C ,  f ile-cfg_f ile,  status-' write'  ) 

CALL  WRITEOUT (10) 
close (10) 
endif 


9999  return 
end 


'  Entry  format  error  --  try  again  ' 

'  *************★**********'**********************' 
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subroutine  trim(char_sti, lengthspec, : 1 , i2 ) 


r  mds  positions  of  the  first  and  last  ncn-blan*  characters 
m  a  string 


integer  length  spec, i 1 , i2, i 
character***)  char  str 
ic  1C  i-1 , lengthspec 

if  (char_st r ( i  :  i )  .  ne . '  ')  then 

goto  15 
er.dif 
c  rnt  i  r.ue 

ic  2  0  i-  i  en.gt  h_spec ,  1 ,  - 1 

if  ( cha r_st  r ( i : i )  . ne  .  '  ')  then 

i2“i 
gttc  25 
endi  f 
cent inue 
ret  urn 
end 


SUBROUTINE  WRITECUT  WRITES  ALL  OF  THE  PARAMETERS  FOR  THE  PROBLEM  IN 
THE  FILE  OPENED  AS  file  -  if  lie  IN  THE  CALLING  PROGRAM. 


THE  N'PMAL I  ZED  PARAMETERS  (a.fc.c....)  ARE  PASSED  IN 
AND  WHEN  WP I TE DOT  IS  CALLED  ARE  ASSUMED  NORMALIZED  T 
WRITE  1CT  RENORMALIZES  TO  WAVELENGTHS  IN  THE  HOST  ME: 
3EF0RE  WRITING  OUT  TO  THE  DATA  FILE. 


HF  COMMON  STATEMENT 
FREE  SPACE  WAVE  DEN 7TH 

IA 


SUBROUTINE  WBITEOUT ( if  lie) 

integer  rviim,  le.  rdim,  iunit 
integer  if :le,  s,  •  n 

rea 1  a , t 

real  f  op, f ops top, inerf  r, unit ,  sig 

real  h ( 2 0 » , eps  r ( 20  I , mu  r (20) , Iosa_tan<20) , mu_t an (20) 
Irgical  dc_gc 1 , do  uc  1 
character*50  ccr.f  ig_f  i  *e 


-r-mcr.  pa  r  am '  md  im,  le 
c  “(uric r. ' geom ■  a ,  fc 

ccmrcr.  f  req  fop.  f  cpstcp,  inerf  r 

ccmrson  layers  /  rdiai,  h,  epsr .  »u  r  ,  1  os s _t  ar. ,  n  _*.  an 
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common/ f iles /conf ig_f i le 
common/ unit s/lunit ,  unit 


write 

(ifile, *) 

'***  input  data  file 

write 

(ifile, *) 

ccnfig  file 

write 

(ifile, *) 

'***  maximum  mode  number 

write 

(ifile, * ) 

mdim 

write 

(ifile, *) 

'***  units  (1-inch,  2-cm) 

write 

( ifile, *) 

lun  it 

write 

(ifile, *) 

'***  x  and  y  dimensicr.3 

write 

(ifile, *) 

a ,  b 

write 

(ifile, * ) 

'***  start, stop  and  incrmt  frequencies  in 

write 

(ifile, *) 

fop, fop st op, incrfr 

write 

(ifile, * ) 

'***  number  of  layers 

write 

(ifile, *) 

rdim 

write 

(ifile, *) 

'***  height,  epsr,  mur,  loss_tan leper) , ' 

,  '  loss_tan (mur )  : ' 

do  i  ■ 

1 , rdim 

wr  i 

te ( if i 1 e. 

* )  h ( i) , epsr (i) , mur ( i ) , loss_tan (i)  ,  mu_t an ( i ) 

enddo 

write 

(ifile, *) 

'***  losses  Delow  and  above  walls 

write 

(ifile, * ) 

do_gc 1 

write 

(ifile, * ) 

dc_-cl 

write 

(ifile, «' 

'  •**  conductivity  of  walls 

write 

< i f i -e, * ) 

sig 

write 

dfile,  *) 

write 

(ifile, *) 

return 

end 


SUBROUTINE  LOGO 

THE  UNIVERSITY  of  MICHIGAN  COLLEGE  OF  ENGINEERING 
RADIATION  LABORATORY 
ANN  ARBOR,  MICHIGAN 


UWVI 

MMMM 
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yv 
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MM 
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MM 

MM 

MM  M 

MM 

vy 

MM  MM 
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MM 

MM 

MM 

MM 

MM 

M  MM 

yy 

M  MM 

MM 

MM 

MM 

MM 

MM 
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MMMMMMMM 

MM 

M  MM 

yy 

MM 

MM 

hW 

MM 

MM 

MM 

MM 

MM  MM 

MM 

MM 

MM 

MMM 

yyyyi 

MMMM 

MMMM 

MMRU 

MMMM 

MMMM 

MMMM 

MMMM 

MMMM 

MMMM 

MMM 

MARCH  11,  1990  VERSION 


WRITE ( * , * ) 
WRITE!*,  •) 
WRITE ( * , * ) 
WRITE <*, •) 
WRITE (*, *> 


//tara/us*rs/vand*v*nt*r/mojr*nt_dir/pfw .  ftn  Pag*  19 


03/12/90  12:13 


WRITE ( * , *) 
WPITE  <*, *) 

k 

WRITE  {*,  *) 
WRITE (*, *) 
WRITE  (*,  ») 
WRITE (  «,  *) 

I 

WRITE  <*, *) 

WRITE  (  *, «) 
WRITE  (  *,  *) 


WRITE  (  * 
WRITE ( * 
4 

WRITE ( * 
i 

WRITE ( * 
WRITE  C 

> 

WRITE ( « 
WRITE ( * 
WRITE ( * 

k 

WRITE ( * 
WRITE (* 
i 

WRITE (* 
WRITE (• 
WRITE ( * 
WPITE ( * 
PAUSE 


NEERING' 


THE  UNIVERSITY  OF  MICHIGAN  COLLEGE  OF  FNG I ' , 

RADIATION  LABORATORY' 

ANN  ARBOR,  MICHIGAN' 


r**********i 


i**********************i 


*k*tk*k*k*k*k*kkk*kkk*k**k*kk*#( 

************************ 


MMMM  MMMM  MMMM  MMMM  MMMM  MMMM  MMMM  MMMM  '  , 

MMMM  MMMM  MMMM ' 


WRITE ( 

* 

t 

•) 

t 

MM  M 

M  MM 

MM 

MM 

MM 

MM 

MM 

MM 

MM 

MM  '  , 

i 

0 

MM 

MM 

MM 

M 

MM' 

WRITE ( 

* 

t 

*) 

t 

MM  M 

M  MM 

MM 

MM 

MMMMMMM 

MM 

MM 

0 

0 

4 

0 

MM 

MM 

MM 

M 

MM' 

WRITE ( 

* 

t 

*) 

0 

MM 

M  MM 

MM 

MM 

MM 

MM 

MM 

MM 

MMMM' , 

& 

0 

MMMMMMMM 

MM 

M 

MM  ' 

WRITE ( 

* 

t 

*) 

0 

MM 

MM 

MM 

MM 

MM 

MM 

MM 

MM 

MM 

MM  '  , 

4 

0 

MM 

MM 

MM 

MMM' 

WRITE ( 

* 

1 

*) 

0 

MMMM 

MMMM 

MMMM 

MMMM 

MMMM 

MMMM 

MMMM 

MMMM  ' , 

MMMM  MMMM  MMMM  MMM' 


************************ 

t**************************************************** 


THIS  PROGRAM  CALCULATES  THE  PROPAGATION  CONSTANT  ', 
OF  THE  HYBRID  MODES' 

EXCITED  IN  PARTIALLY-FILLED  WAVEGUIDES' 


T.  EMI LIE  VAN  DEVENTER  —  AND 
KATEHI' 


MARCH  11,  1990  VERSION' 


LINDA  P.  B. '  , 


RETURN 

END 


ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff 

functions 

ffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff 
complex  function  eqn_c(kz) 

complex  sum, f sum, kz 
common /t rick/ f sum 


sum  -  cmplx(0.,0.) 
call  hybrid (kz, sum) 
eqn_c  ■  sum 
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f 3um*sum 


return 

end 

cffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffffff 
complex  function  twg (arg) 
complex  arg,  j 


■i  «  <~mr>  1  x  (  n  n  1.0' 
if (aimag (arg)  .gt.  45.0)  then 
twg  -  j 

else  if  (aimag (arg!  .It.  -45.0)  then 
twg  =  -j 
else 

twg  -  csin(arg)  /  ccos(arg) 
endi  f 
return 
end 

cf f ff ffff ffffff ff fffff f ffffffffffffff ffffffffffffff f ff ffff ff fff ff fff fff 
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